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FREQUENCY OF ENCOUNTER OF AIRCRAFT 
IN A RANDOM HORIZONTAL FIELD 

John D. Bird and Kathryn A. Smith 
Langley Research Center 

SUMMARY 

Calculations have been made of the frequency of encounter as a function of azimuth 
of encounter of a passing aircraft with the aircraft in a random planar horizontal field. 
The field aircraft were all moving at a constant speed but in random directions. These 
calculations included the total frequency of encounter with the aircraft of the field and 
the frequency of encounter with those aircraft of the field which were encountered in the 
fore quadrant, in the lateral quadrants, and in the rear quadrant; the calculations were 
made for various speed ratios of the field aircraft and the passing aircraft. In the 
cases studied, the passing aircraft was usually approached in its fore quadrant. Rela- 
tively few of the aircraft of the random field were approached from the rear because a 
speed subtraction effect reduced the encounter velocity in the overtaking situation. The 
frequency of encounter with the field aircraft which were approached from the lateral 
quadrant varied little with the speed ratio of the aircraft because in this case there was 
no subtractive or additive effect of the velocities. 

INTRODUCTION 

The study of chance encounter of aircraft is important to both civil and military 
aviation because of collision hazard and is important also to military tactical combat 
situations. Attention has recently been given to collision probability and to collision- 
avoidance techniques in air traffic control studies. (See refs. 1 and 2.) Terminal-area 
collision hazard has been the matter of greatest concern in this work, and both experi- 
mentally determined and statistically modeled collision probabilities have been studied. 

This investigation provides calculations of the number of aircraft encountered by 
one passing aircraft traveling at a constant speed and in a constant direction through a 
uniform and unlimited horizontal planar field of aircraft moving randomly in direction at 
a second constant speed. The number of encounters per second per radian of the passing 
aircraft with the aircraft of the field was calculated as a function of the azimuth of 
encounter relative to the direction of flight of the passing aircraft. No form of air traffic 
control was considered. The integral of this quantity around the aircraft (the total 



frequency of encounter) was also determined. The frequency of encounter of field air- 
craft by the passing aircraft is listed in categories corresponding to (a) those aircraft 
of the field approached from the fore quadrant, (b) those approached from the lateral 
quadrants, (c) those approached from the aft quadrant, and (d) the total aircraft encoun- 
tered. (See fig. 1.) Calculations to determine the effect of relative velocities were 
made, and the results were normalized with respect to the number density of the aircraft 
in the random field. Considerations of rarified gas dynamics were used in these calcu- 
lations, but a distribution function which is random in direction and uniform in magnitude 
of velocity was employed rather than the Maxwellian distribution common to gas dynamics. 
This distribution function seems more appropriate than the gas-dynamics distribution 
function for illustrating the significance of the various parameters of the aircraft- 
encounter problem. The principal features of these results are discussed in regard to 
collision hazard and aerial combat tactics. (See ref. 3.) 

SYMBOLS 

aircraft passing at constant speed and direction through field of aircraft 
moving randomly in direction at another constant speed 

random field of aircraft 

constant 

number of aircraft of field B that enter circle of encounter around 
aircraft A 

frequency at which aircraft cross circle of encounter 

frequency at which aircraft cross circle of encounter as function of 6 
radius of encounter around aircraft A 
element of circle of encounter of aircraft A 
vector form of ds 

speed of aircraft A moving along axis Y relative to axes x' and y' 
speed of aircraft of field B relative to axes x' and y' 


A 

B 

k 

N 

dN 

W 

d 2 N 
dt d0 

r 

ds 
ds 
u 
v T 
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V 


vector velocity of aircraft in field B relative to ds 


X,Y 

x,y 

x’,y’ 

8 

6 o 

e 

9 

<t> 

h 

0 * 


speed of aircraft of field B with respect to ds 
axes fixed to aircraft A with Y alined with inertial axis y* 
axes fixed to aircraft A at angle Q to axes X,Y 
inertial system of axes 

direction of flight of aircraft of field B measured from axis x 

specific value of 5 

angle measured around aircraft A 

density of aircraft in field B 

direction of flight of one aircraft of field B measured from axis x* 
lower limit of (p 
upper limit of cp 

specific values of cp relating to integration limits 
defined value of <p 

THEORETICAL DEVELOPMENT 


Consider an aircraft A which is moving in a given direction at a constant speed u 
through a uniform and unlimited horizontal planar field of aircraft B which are moving 
randomly in direction but at a constant speed v' . (See fig. 1.) Calculation of the fre- 
quency at which the aircraft of field B reach a circle of encounter of radius r encom- 
passing and centered on aircraft A under various conditions is now desired. Specifi- 
cally, calculations of the frequency at which the aircraft of field B reach r at various 
angles about aircraft A, the total frequency at which aircraft of field B are encountered 
by aircraft A, and the frequency at which aircraft of field B are encountered from the 
fore, lateral, and rear quadrants at various angles about aircraft A are desired. The 
aircraft of field B were considered as individual points in the field. 
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The x- and y-axis system of figure 2 is fixed relative to aircraft A for a given 8 
and r; the x- and y-axis system moves along the direction Y with speed u relative to 
the inertial system x',y\ The aircraft B move at speed v' in a direction 0 rela- 
tive to x'.y' and have a field density d 0 where this field density is the number of 

2 7 r 

aircraft B per unit of area in angle d 0. 

The velocity v of the encountered aircraft of field B relative to the element ds, 
the vector element ds, and the field density d0 determine the frequency at which 

" 7T 

aircraft enter the encounter circle. (See fig. 3.) This frequency is 


d far) - If A(p ds ' v 

or the rate per unit of 8 is 


d 2 N _ pr 
dt d6 2 tt x 


d 0 


where 


v = v' cos (8 - 0) - u sin 8 

An integration is set up over the angle 0 for the limits to be considered for all 
the aircraft of field B which approach the encounter circle from the outside. With 
certain changes, this formulation gives 


d 2 N 
dt d e 




\ 


^ cos (e - 0 ) - 


sin 9 


d 0 


where 


h = 8 - 6 0 
0 U = e + 6 0 

The quantity ±6 Q defines an angle enclosing the aircraft of field B that are under 
consideration for possible encounter with aircraft A. (See fig. 2.) An angle 6 0 of 180° 
encompasses all the aircraft of field B; an angle 6 Q of 45° encompasses only those 
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aircraft of field B which are approached from the rear. When the integrand is negative, 
k = 1; when the integrand is positive, k = 0. This condition insures that only those air- 
craft in the field B which approach the encounter circle from the outside are included. 
The completed integration gives 


d 2 N 
dt d0 



0 sin (p - sin 9 cos (p) - <p sin 6 


0u 


In evaluating this expression, the conditions where the integrand v X/ /u is positive 
must be determined to eliminate the unwanted count of aircraft in field B. The following 
expression is examined: 


= -5- cos (e - 0) - sin 6 


where 



(0 - 6 0 ) § <p § (e + e 0 ) 


0 g 6 0 g 7T 


Thus, v /u is seen to be symmetrical about cp= 6 with a maximum value at <p = 9, 

^ V v v f 

at which point -jp = — - sin 9. At the lower and upper limits of cp, 
v x v' 

— = — cos 6 0 - sin 0 

There are three cases relating to the conditions for which v x/ /u is positive or 
negative. (See fig. 4.) These cases are: 

(1) When ^ — sin 0> 0 and cos 6 n - sin 0 > 0, v /u is positive for all (p 
from to cp n and k = 0 in the expression for d^N/dt d0. In this case 

d 2 N/dt d0 = 0. 

TT t « T 

(2) When — - sin 0 § 0 and — cos 6 Q - sin 0 s 0, v x/ u changes sign twice 
between <p^ and cp u and is negative at <p ^ and cp^. The points where v/u changes 
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sign are given by equating v x/ /u to zero and solving for <p^ and <^. This process 
gives 


cp 1 = e - 

(j>2 = d + cf>* 


where 


(p* = cos~l^- sin 6j 

In this case the expression for d 2 N/ dt d0 is evaluated from <p ^ to <p^ and <pn 
to <p u with k = 1. The constant k is zero from to (p^. 

-*r f y t 

(3) When - sin 6 < 0 and — cos 6 n - sin 6 < 0, v /u is negative for 

all (p from (p ^ to <P u and the expression for d^N/dt d0 is evaluated with k = 1 
from to to obtain the result. 

The limits of the function sin 6 for these three cases are shown in figure 5, and 
the evaluation of the quantity d 2 N/dt d0 for each of these cases is given here. 

Case (1): 


-dfN =0 

dt d 6 %.;■ 

Case (2): “ 

Sll = 2 (fr)[/r( sl " 5 ° " sln **) ' ( 5 ° ' **) sln 9 

Case (3): 


d 2 N 
dt d6 




sin 6 q - 


6 0 sin 6\ 


The total rate of encounter of aircraft A with the aircraft of field B dN/dt is 
obtained by integrating d 2 N/dt dd over the angle 6 completely around aircraft A 
with a numerical quadrature routine. In this evaluation of dN/dt, various possibilities 
which relate to the exclusion of situations where sin 6 exceeds its normal limits of ±1 
must be considered. These situations are given in table I. 
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CALCULATIONS 


Calculations of the encounter per second by an aircraft A moving at a constant 
speed in a field of aircraft B, all of which are moving in random directions but at a 
constant speed, were made for a field density factor of pru/27T = 1 per second for speed 
ratios of the aircraft B to aircraft A of 0, 0.5, 0.8, 1.0, 1.2, and 2.0. The frequency 
of encounter of aircraft A with the aircraft of field B, in units of encounters per second 
per radian, was calculated as a function of the azimuth position of encounter relative to 
the direction of flight of aircraft A. Calculations included the encounters in which air- 
craft A approached the aircraft of field B from the fore quadrant, from the lateral 
quadrants, and from the aft quadrant. Calculations were also made for the total frequency 
of encounter in each of these cases, that is, the integral of the encounters per second over 
the azimuth angle 6 of encounter. 

The calculations of frequency of encounter in encounters per second per radian are 
plotted as a function of azimuth angle 6 of encounter in figures 6 to 9 for various speed 
ratios. The total frequency of encounter is listed in table II as a function of v'/u and 
is plotted in figure 10. 


RESULTS 

Figures 6 to 9 show the frequency of encounter d^N/dt d 8 of the aircraft in a 
random planar field B with a passing aircraft A as a function of the azimuth of 
encounter 6 with aircraft A for various speed ratios v’/u for four cases. These 
cases correspond to the total frequency of encounter and to the frequency of encounter of 
those aircraft of field B which are encountered in the fore quadrant, lateral quadrant, 
and aft quadrant. A general pattern of increasing frequency of encounter with increasing 
values of v'/u can readily be seen in all the figures. This trend indicates an increased 
frequency of traverse of a given area in the field as the speed v’ is increased. Most 
of the encounters are made in the fore quadrant of aircraft A; this situation particularly 
applies in the case of the aircraft of field B which are encountered in the lateral quad- 
rants (fig. 8). Relatively few aircraft of the field B are encountered in the aft quadrant 
because of the slow rate at which the passing aircraft A overtakes the aircraft of the 
field B (fig. 9). 

The functions of figures 6 to 9 were integrated to give the total number of encounters 
experienced by the passing aircraft A for the various cases. These results are plotted 
against speed ratio v’/u in figure 10 and are listed in table II. The limiting case where 
v'/u goes to zero is also included in figures 6 to 10 and table II. The character of the 
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curve of figure 10 for those aircraft of field B which are encountered from the lateral 
quadrants is somewhat different from the curves for those aircraft approached from the 
fore or aft quadrant. At a speed ratio v'/u of about 1.2, the curve for those aircraft 
of the field B which are encountered from the rear goes to zero and remains zero for 
larger values of v'/u. The curve for those aircraft which are encountered from the 
sides is less affected by speed ratio than the other curves. The curve for the aircraft 
which are approached from the front increases rapidly as the ratio v'/u is increased. 
This effect results from a more rapid sweep of the encounter area by the aircraft of the 
field as their speed increases. The results presented here may be employed to derive 
frequency of encounter for specific situations by calculating pru/27r and then multiply- 
ing by the encounter frequency for the appropriate value of v’/u. 

CONCLUDING REMARKS 

Calculations were made of the frequency of encounter as a function of azimuth of 
encounter of a passing aircraft with the aircraft of a random planar horizontal field. 

All field aircraft were moving at a constant speed but in random directions. These calcu- 
lations included the determination of the total frequency of encounter by the passing air- 
craft with the aircraft of the field and the frequency of encounter by the passing aircraft 
with those aircraft of the field which were encountered in the fore quadrant, the lateral 
quadrants, and the aft quadrant for various speed ratios of the field of aircraft and the 
passing aircraft. The results of these calculations are presented in terms of the frequency 
of encounter as a function of the azimuth of encounter around the passing aircraft for the 
various speed ratios; the results are also presented as integrated or total frequencies 
of encounter for the various conditions. 

The encounters with the passing aircraft were largely in the fore quadrant of the 
passing aircraft; relatively few of the aircraft of the field were approached from the rear 
as a result of a speed subtraction effect which reduced the encounter velocity in the over- 
taking situation. The frequency of encounter with the aircraft in the field which were 
approached from the sides varied little with the speed ratio of the aircraft because there 
was no subtractive or additive effect of the velocities. Basically, the largest number of 
encounter frequencies tended to occur in the fore quadrant both of the passing aircraft and 
of the aircraft of the random field. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
January 14, 1976 
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TABLE I.- POSSIBILITIES TO BE CONSIDERED IN EVALUATION OF 
TOTAL FREQUENCY OF ENCOUNTER dN/dt 


Condition 

7T C0S 8 Q S-1 

-1 < TT cos 6 o < 1 

V T 

— cos 6 > 1 

u o - 


Not possible 

All cases are 
possible 

Not possible 

Hsl 

U “ 

Only case (2) 
is possible 

Only case (1) and 
case (2) are 
possible 

Only case (1) 
is possible 


TABLE II.- FREQUENCY OF ENCOUNTER BY PASSING AIRCRAFT 
WITH AIRCRAFT OF UNIFORM PLANAR FIELD 
|pru/27r = l.oj 


Ratio of speed of field 
to speed of passing 
aircraft, v'/u 

Encounters per second, dN/dt 

Approached 
from fore 
quadrant 

Approached 
from lateral 
quadrants 

Approached 
from aft 
quadrant 

Total 

0 

3.14 

6.28 

3.14 

12.56 

.5 

5.69 

6.43 

1.25 

13.37 

.8 

7.57 

6.66 

.45 

14.68 

1.0 

9.00 

6.87 

.12 

15.99 

1.2 

10.68 

7.14 

.02 

17.84 

2.0 

17.77 

8.96 

0 

26.73 
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-sin 0> 0 


cos 6 


-sin 0> 0 


— -sin 050 
u 

— cos 6 -sin 0^0 
u o 


— — sin 0 < 0 
v' 

— cos 6 -sin 0 < 0 
u o 


Figure 4.- Three cases that relate to sign of v /u. 




Case (1) : sin 0 < — cos 6 
u o 

Case (2) cos 6 ^ sin 0£— 
u o u 

Case (3) : sin 0>-^ 

Figure 5.- Limits of sin 9 for three cases relating to v x j u. 
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Y 
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Figure 7.- Frequency of frontal encounter with aircraft of field B 
of encounter with passing aircraft for various speed ratios, 


as function of azimuth 
pru/27r = 1.0. 
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dN 

dt * 

N encounters 
per sec 



Figure 10.- Frequency of encounter by passing aircraft with aircraft of uniform 

planar field, pru/2ir = 1.0. 
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